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Abstract 

The three-dimensional bimodal random-field Ising model is studied via a new finite temperature 
numerical approach. The methods of Wang-Landau sampling and broad histogram are imple- 
mented in a unified algorithm by using the N-fold version of the Wang-Landau algorithm. The 
simulations are performed in dominant energy subspaces, determined by the recently developed 
critical minimum energy subspace technique. The random-fields are obtained from a bimodal dis- 
tribution, that is we consider the discrete (±A) case and the model is studied on cubic lattices with 
sizes 4 < L < 20. In order to extract information for the relevant probability distributions of the 
specific heat and susceptibility peaks, large samples of random-field realizations are generated. The 
general aspects of the model's scaling behavior are discussed and the process of averaging finite-size 
anomalies in random systems is re-examined under the prism of the lack of self-averaging of the 
specific heat and susceptibility of the model. 
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I. INTRODUCTION 

The random-field Ising model (RFIM) [ill is one of the most studied glassy magnetic 

n n n 

models [2, 3, 4]. The 3D RFIM consists of Ising spins Si on a simple cubic lattice, governed 
by the Hamiltonian: 

<i,j> i 

where J > is the interaction constant and hi are quenched random-fields, obtained from a 
bimodal distribution P(hi) = \5{hi — A) + ^8(hi + A). A denotes the disorder strength, also 
called randomness of the model. Although nowadays it is believed that the phase transition 
from the ordered to the disordered phase of the model is of second-order, a complete set of 
critical exponents fulfilling a widely accepted set of scaling relations has not been established. 
In fact, there is a strong disagreement in literature concerning the overall thermal and 
magnetic behavior of the model 0, Q|. This may be due to a mistaken comprehension of 
some theoretical concepts in random systems, such as the concept of averaging that will be 
discussed below. 

The rest of the paper is laid out as follows. In Section [Til we present the numerical 
schemes utilized for the study of the RFIM. The process of averaging finite-size anomalies 
in random systems and the significance and implications of the non trivial property of the 
lack of self-averaging of the specific heat and susceptibility of the model, are discussed in 
Section ITTT1 Finally, we summarize in Section ITVl 

II. NUMERICAL TECHNIQUES 

Numerically the RFIM has been approached using traditional [?], Q] but also more sophis- 
ticated Monte Carlo techniques Q|. However, the nature of the model demands enormous 
computer resources. Furthermore, in order to get a good estimate of the mean properties 
of the system, it is necessary to repeat the simulations for a large number of realizations 
of the random-fields. Here, the numerical procedure concentrates on the determination of 
the density of states (DOS) G(E) of the model and on the corresponding thermodynamic 
quantities. 

For the application of the Wang-Landau (WL) algorithm jlO] in a multi-range approach 
we follow the N-fold description of Schulz et al. [11]. The random walk is not allowed to 



2 



move outside of any particular subrange, and we always increment the histogram H(E) — > 
H(E) + 1 and the DOS G(E) —>■ G(E) * fj after a spin-flip trial. Here, fj is the value of 
the WL modification factor / [10] at the jth iteration, in the process (/ — > f 1 ^ 2 ) of reducing 
its value to 1, where the detailed balance condition is satisfied. In all our simulations the 
control parameter takes the initial value: fj = % = e ~ 2.71828..., while when starting a new 
iteration it is changed according to the sequence fj + i = ^/fj,j = 1, 2, 20 [h], Q|- For the 
histogram flatness criterion we use a flatness level of 0.05. The accumulation of numerical 
data for the app.ication of the broad histogram (BH) method of Oliveira et al Q and a.so 
the updating of appropriate (E, M) histograms is carried out in the final stage of the WL 
process, by using the N-fold iteration j = 12 — 20 jl^ . The approximation of the DOS, in the 
last WL iteration, Gwl{E), and the high-level (j > 1) WL (E, M) histograms, H WL (E, M), 
are then used to estimate the magnetic properties in a temperature range, which is covered, 
by the restricted energy subspace (Ei,E 2 ) as: 

(Mn) E E (M n ) E G(E)e-^ Y, EW2) {M n ) E , WL G WL {E)e-^ 

{ } ~ E E G(E)e-^ ^ E E ^E u E 2) GwL(E)e-^ {) 

The micro canonical averages (M u )e are obtained from the Hwl{E,M) histograms as: 

(m") e = {m^ )e , wl ^y.^ Bw „ l{e :^ ) 

M H WL{E) 

H wl (E) = Y j H W l(E,M) (3) 

M 

and the summation in M runs over all values generated in the restricted energy subspace 
[Ei, E<z) [14]. Similarly we obtain the microcanonical estimators necessary for the application 
of the BH method, using the well-known broad histogram equation [ljj: 

G(E)(N(E, E + AE n )) E = G(E + AE n )(N(E + AE n , E)) E+AEn (4) 

where N(E, E + AE n ) is the number of possible spin flip moves from a microstate of energy 
E to a microstate with energy E + AE n , which are known during the N-fold process. 

For a particular random-field realization the specific heat and its peak are easily ob- 
tained with the help of the usual statistical sums. The critical minimum energy subspace 
(CrMES) scheme [12I . Q] uses only a small but dominant part {E-, E + ) of the total energy 
space (E m i n , E max ) to determine the specific heat peaks. Let E denotes the value of energy 
producing the maximum term in the partition function at the pseudocritical temperature 



(corresponding to the specific heat peak) and S(E) = In G(E) the microcanonical entropy. 
Then the CrMES approximation is defined by the following equations: 

C L (E^E + ) = N~ l T- 2 \z- 1 Y,E 2 ex V [^(E)]- ( Z' 1 £ Eexp [5(E)] ] 1 (5) 



E+ 

$(E) = [S(E) - PE) - [S(E) -{3E\, Z = J^ ex P ( 6 ) 
where (E_, E + ) is the minimum dominant subrange, satisfying the following accuracy crite- 



rion: 



C L (E_,E + ) 1 
C 1 L^E m i n , E max ^j 



< r (7) 



with r = 1CT 6 . Note that, the above accuracy is extremely demanding compared to the 
statistical errors produced by the DOS method (i.e. the WL method) and to the large 
sample-to-sample fluctuations of the RFIM that will be discussed below in Section IIHI 

Using an ensemble of macroscopic samples of size L corresponding to different random- 
field realizations we have applied the described scheme in a broad energy(magnetization) 
space (total CrME(M)S of the ensemble) that covers the overlap of the dominant en- 
ergy(magnetization) subspaces for all realizations of the ensemble. This practice has the 
advantage that the approximation of the specific heat and susceptibility for a particular 
random-field is accurate in a wide temperature range, including its pseudocritical tempera- 
ture. Despite the strong fluctuations of the energy value corresponding to the maximum term 
of the partition function Z, the union of the CrME(M)S for large samples of random- fields 
is, in any case, a quite small subspace. 



III. AVERAGING FINITE-SIZE ANOMALIES. LACK OF SELF- AVERAGING 

For a disordered system one has to perform two distinct kinds of averaging. Firstly, for 
each random-field realization the usual thermal average has to be carried out and secondly 
one must average over the distribution of the random parameters. The latter makes it clear 
that large ensembles of random-fields must be generated in order to estimate properly the 
mean properties of the system. Following the methods described above in Section [III the 
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thermal average for the specific heat is given by Eq. (j3j), while the susceptibility x reads as: 

X = ^{(M 2 )-(M) 2 } (8) 

with N = L 3 . Let C m (T) and Xm(T) denote the specific heat and susceptibility of a par- 
ticular random-field realization m in an ensemble of M realizations (m = 1, 2, M). The 
corresponding pseudocritical temperatures T£(C m (T)) and T£(% m (T)) depend on the par- 
ticular realization of the random-field and for large values of the randomness A, they are 
strongly fluctuating quantities. The locations of the specific heat and susceptibility peaks 
may be then denoted bv (Cm, T L,c-,m) and (Xm.^m): respectively 

In previous studies |7|, [8], the averaging process over a large number of random-fields has 
been carried out on the averaged curve of the specific heat or susceptibility without first 
raising the question of whether this averaged curve is the proper statistical representative of 
the system. Specifically, the following sample averages have been considered for the specific 
heat and susceptibility [3, Q] : 

1 m 1 M 

[ c u = ^J2 c m(n fcu^ExmCn ( 9 ) 

m=l m=l 

The finite-size scaling behavior of the peak of these averaged curves was then studied by 
assuming that the maxima [C]* av = ([C] av )* and [x]* av = (\x\av)* obey power laws (for details 
see Refs. [l, sj]). Note that, these averaged curves ([C] „ and [x]o«) are very sensitive to the 
property of self-averaging due to the fact that the corresponding thermodynamic quantities 
are characterized by broad distributions in the thermodynamic limit. 

In this work, in addition to the above averaging expressions, we study the sample-averages 
of the individual specific heat and susceptibility maxima, defined by: 

j M I M 

m=l m=l 

These mean values, together with the corresponding peaks of the averaged curves of Eq. (Q, 
are shown in Figs. [Tf2l In our simulations we used an ensemble of M — 1000 random-field 
realizations for L < 12 and M = 500 for L = 14 — 20. To quantify the sample-to-sample 
fluctuations of the specific heat (susceptibility) peaks we define the standard deviation of 
C* m (Xm)> over a sam ple °f M random-field realizations as a(C^) (&(Xm))- This important 
parameter will be illustrated in our figures as error bars, but should not be in any case 
confused with the existing statistical errors. 
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FIG. 1: Finite-size behavior of the averages [C^] a „ and [C]* v for the case A = 2, for both the WL 
and BH methods used. The error bars represent the sample-to-sample fluctuations (see text). The 
behavior of [C]* av appears as a random fluctuation around the value 0.815, as shown by the dotted 
line. 

Fig. [T] illustrates the finite-size behavior of the peaks of the average [ClJ a .„ and that of 
the averaged curve [C]* v , defined above for the case A = 2 for the two methods employed, 
i.e. the WL and BH methods. As discussed above, the error bars quantify the large sample- 
to-sample fluctuations of the specific heat peaks for the two methods employed (note that 
the error bars with the larger cap- width always refer to the WL method). From Fig. [I] it 
is apparent that, while the sample mean averages [ClJ a „ admits of finite- size scaling, the 
behavior of [C]* v appears as a random fluctuation around the value [C]* v ~ 0.815, as shown 
by the dotted line in this figure. In analogy with Fig. [U Fig. [2] shows the corresponding 
susceptibility quantities, also for the case A = 2. Again, the sample-to-sample fluctuations 
are very large and the behavior of [x\* av is clearly distinct from that of [XmW Inspecting 
Figs. [T1I21 on a comparative basis we observe that the deviations between the WL and BH 
methods are more pronounced in the thermal case (Fig. [T]) and their difference is noticeable 
for L = 20. This difference represents the order of the statistical errors of our scheme. 
Although these errors are still small compared to the sample-to-sample fluctuations, they 
raise doubts whether the total number of WL iterations {j final = 20) is sufficient for the 
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study of large lattice sizes. 
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FIG. 2: Finite-size behavior of the averages [Xm]<™ an( l \x\av f° r the WL and BH methods, also for 
the case A = 2. 



From the discussion above, it is clear that when studying random systems the only 
meaningful objects for investigating the finite-size scaling behavior are the distributions 
of various properties in ensembles of several realizations of the randomness. Hence, it is 
important to be able to ascertain to what extent are the results obtained from an ensemble 
of random realizations representative of the general class to which the system belongs. The 
answer hinges on the important issue of self-averaging. If a quantity is not self-averaging, 
we talk about lack of self-averaging and the process of increasing L does not improve the 
statistics. In other words, the sample-to-sample fluctuations remain large. The problem of 
self- averaging in the 3D RFIM has been a matter of investigation over the last years [15|. A 
common measure characterizing the self-averaging property of a system based on the theory 
of finite-size scaling has been discussed by Binder [16| and has been used for the study of 



some random systems 



13, 



width quantity, defined as: 



181 ] . This measure inspects the behavior of a normalized square 
Vn 



Rc 



[QY 



(ii) 



where Vq = [Q 2 ] — [Q] 2 is the sample-to-sample variance of the average [Q]. Here, Q is used in 



respect of the specific heat C and the susceptibility \- According to the literature 
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when the ratio Rq tends to a constant, the system is said to be non self- averaging and the 
corresponding distribution (say P(Q)) does not become sharp in the thermodynamic limit. 




FIG. 3: Temperature variation of the ratio i2|gj defined in the text, for both the specific heat 
and susceptibility. 

Using our notation we may define the ratio Rq for the specific heat and susceptibility in 
two explicit forms, one for the case of the averaged curve [Q] av - 

[c] - ~ (icur v ~ {[x)c ' 2 



lav i 



and one for the case of the average [Q 



m)av 



D V lC* m ]av _ ( °{C*m) \ 2 p t V lX*m]av _ ( °{Xm) V ( 1 o\ 

\\ym\av) \y^m\av / \[Xm}av) \[Am\av / 

In Fig. [3] we present the behavior of the ratio R'\Q\ av = R[Q] av / 'R*Q] av i where R*q^ av = 
max{i?[Q] at) } as a function of the temperature T, for L = 16 and A = 2. The solid line 
corresponds to the specific heat (Rip,^) while the dotted line to the susceptibility C^Li )• 
In this figure only the results of the WL method are presented, since our intention was to 
identify the temperature variation of the non self-averaging property of the averaged specific 
heat and susceptibility defined in Eq. (Q. Indeed, from Fig. [3] we observe that for tempera- 
tures close to the critical, the ratio R[Q] av is maximized indicating strongly non self-averaging 
behavior for both quantities. In Fig. H]we consider the behavior of the ratio R[Q* m ] av of the 
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FIG. 4: Illustration of the ratio R[q* ] av of the specific heat and susceptibility, for both the WL 
and BH methods, for randomness A = 2. Clear saturation to a limiting non-zero constant value: 
R ]c * i -» 0.28 and R w * i -» 0.57. 

specific heat (R[c^] a v) an d susceptibility (-R[ X £j a J as a function of the linear size L, for both 
the WL and BH methods. Both ratios seem to tend to a constant non-zero value, namely 
R[c^\av ~^ 0-28 and R[ x * m ] av — > 0.57, confirming the above mentioned lack of self-averaging 
of the specific heat and susceptibility of the model. 

IV. SUMMARY AND OUTLOOK 

The numerical strategy applied in this paper enabled us to perform extensive finite- 
temperature simulations and extract valuable information for the generic behavior of the 
RFIM. It was shown that the various definitions of the apparent finite-size anomalies may 
not be equivalent. Our analysis revealed that the behavior of the mean [Qmlav is clearly 
distinct from that of [Q]* v and that this is directly connected to the lack of self-averaging 
of the model. More work needs to be done towards this direction, so that the subtle matter 
of self-averaging in the RFIM is fully clarified and understood. One of our future plans, 
is the verification of the above results by studying the model for larger lattice sizes and 
various randomness values. In any case, the present study puts forward some new ideas and 
an efficient unified implementation of the DOS methods, suitable for the study of random 
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systems. 
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